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Abstract. Knowing the symmetries of a polyhedron P can be very useful for 
practical polyhedral computations as well for analysis of the structure of P. 

We study the groups preserving the linear, projective and combinatorial struc- 
ture of P. In each case we give algorithms to compute the symmetry group and 
discuss some practical experiences with these algorithms. 

Our focus here is on R"; we observe that some of the central notions do 
not admit a straightforward generalization to point configurations in complex 
projective spaces. 
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1. Introduction 

Symmetric polyhedra occur frequently in diverse contexts of mathematics. Poly- 
hedra are central to the theory of Mathematical Optimization (Mathematical Pro- 
gramming), and main objects of study in linear and integer linear programming. 
Frequently studied symmetric polyhedra in applications such as transportation lo- 
gistics or machine scheduling have names like "Travelling Salesman", "Assign- 
ment" or "Matching". For these and further examples we refer to [32] and the 
numerous references therein. Polyhedra play also prominent roles in other parts 
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of mathematics, e.g. in algebraic geometry, and in particular in the theory of toric 
varieties [6]. 

For the analysis of high dimensional polyhedra it is important to know their 
symmetries. Furthermore, for many important tasks in polyhedral computations, 
such as linear and integer linear programming, the representation conversion prob- 
lem or volume computations, symmetry exploiting techniques are available (see 
[24] and [5]). Even commercial optimization software like [41] and [42] include 
some techniques for symmetry exploitation by now. To a large extent, the used 
methods depend on the kind of symmetry that is available and how it is presented. 
For instance, if we know the affine symmetry group of a polyhedron from a linear 
programming problem, we can reduce the problem dimension if the utility function 
is invariant under the group. In contrast to the full combinatorial symmetry group, 
the affine symmetry group has the advantage that it can practically be computed 
using only partial information, for instance, from a description of the polyhedron 
by linear inequalities. 

In this paper we provide an overview on how to determine different symmetries 
of a polyhedron computationally. We provide a collection of computational recipes 
for the main polyhedral symmetry groups of interest. Our general philosophy is the 
translation into a problem of determining all the combinatorial automorphisms of 
a colored graph. Although these graph automorphism problems are not completely 
understood from a complexity theoretical point of view (see [16]), there exist so- 
phisticated software tools for their practical solution [33, 36, 37]. All algorithms 
explained here are available within the GAP package polyhedral [34]. The lin- 
ear group of a polyhedron can also be computed with the C++ tool SymPol [40]. 
Similar computational tasks for special classes of lattice polytopes are performed 
by PALP [38, 19]. PALP arose from the needs of the project to classify certain 
lattice polytopes called reflexive polyhedra [17, 18]. 

The paper is organized as follows. In Section 2 we define the three most impor- 
tant polyhedral symmetry groups and explain some of their relations among each 
other: The linear symmetry group, the projective symmetry group and the com- 
binatorial symmetry group. In the following sections we consider each of them 
separately. We start with the linear symmetry group in Section 3. In Section 3.2 
and Section 3.3 we deal with particular subgroups of the linear symmetry group 
that are important to Integer Linear Programming and applications in Computa- 
tional Geometry of Positive Definite Quadratic Forms. In Section 4 we consider 
the projective symmetry group. We give a new characterization and from it derive 
a new algorithm for its computation. In our last Section 5, we consider a practical 
approach to the computation of the combinatorial symmetry group, which contains 
all the other groups. 
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2. Polyhedral Symmetry Groups 

A polyhedral cone C in the vector space M" is defined as the set of vectors 
satisfying a finite number of linear (homogeneous) inequalities. By the Farkas- 
Minkowski-Weyl Theorem there exists a second (dual) description 

C = {Aif 1 + . . . + XpVp : Aj G M>o} 

with generating vectors vi, . . . ,Vp and extreme rays 

Ri = M>QVi. 

Whereas the extreme rays of C are uniquely determined, the generating vectors are 
not. 

A face of a polyhedral cone C is the intersection of C with a supporting hyper- 
plane, that is, with a hyperplane {x G M" : a*x = O} for some a G M" such that 
C is contained in the halfspace {x G M" : a^x > O}. Faces are partially ordered 
by setwise inclusion, which gives a combinatorial structure that is called the face 
lattice of C. The dimension of a face is defined as the dimension of the smallest 
linear subspace containing it. Without loss of generality, we assume that C is full- 
dimensional, i.e. has dimension n and that it does not contain a non-trivial vector 
space. Facets of C are faces of dimension n — 1. The set of facets and the set of 
extreme rays are uniquely determined by C and the problem of passing from one 
description to the other is called the dual description problem. 

A polytope P is the convex hull of a finite set of vectors (wj)i<j<j\/ in M". 
By considering the vectors f- = and the polyhedral cone C defined by 

(t>-)i<j<M, we can actually embed P into C and translate the notions introduced for 
polyhedral cones to polytopes. For more information and background on polyhedra 
and polytopes we refer to [26]. 

By we denote the set of fe-dimensional faces (A;-faces) of C. Such /c-faces 
are identified with the set of extreme rays contained in them. The combinatorial 
symmetry group Comb(C) of C is the group of all permutations of extreme rays 
that preserve Tk for all < A; < n — 1. So in particular, Comb(C) is a subgroup 
of the symmetric group Sym(p) on p elements. 

It is well known that Comb(C) is actually determined by Tn-\- Comb(C) is 
isomorphic to the automorphism group of the bipartite facet-ray-incidence graph. 
Indeed, there is generally no simpler way to compute Comb(C), as this problem 
is graph isomorphism complete even for simple or simplicial polytopes (see [14]). 
For the construction of this bipartite graph we must know both, the extreme rays 
and the facets of C. However, in practice, usually only one of these descriptions is 
known. 

Let GL(C) be the group of invertible matrices A G GL„(M), which induce a 
projective symmetry of C, that is AC = C. The action of GL(C) on C induces 
a permutation of the extreme rays of C. We call the resulting permutation group 
the projective symmetry group Proj(C) of C, which is a subgroup of Comb(C). In 
other words, elements of Proj(C) are permutations a G Sym(p) for which there 
exist matrices A G GL(C) such that ARi = Rfj(j^ for 1 < i < p. 
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^'i P2 P3 Pa 



Figure 1. 4 polytopes in Euclidean space. The classes of 
equivalence under linear, projective and combinatorial equiv- 
alence are ({Pi, P2}, {P3}, {A}), ({Pi, P2, P3}, {A}) and 
({Pi, P2, P3, Pa}) respectively. 



Note, however, that GL(C) and Proj(C) are not isomorphic since the kernel K 
of the homomorphism GL(C) — )• Proj(C) is non-trivial. It contains for instance all 
dilations. In group-theoretic terms, GL(C) = K x Proj(C), cf. Theorem 4. In a 
more general setting of symmetries of a configuration of p points Ri in a projective 
space over a field (e.g. over C), the analog of the group Proj(C) does not need 
to exist, as K does not need to be split from GL(C). See Remark 5 for a simple 
example of the latter, and [30, 13] for the related group-theoretic notions. 

Fixing a set u = {uj}i<i<p of generators for the extreme rays {Ri}i<i<p, we 
obtain in a similar way the linear symmetry group Lin^(C) of C (with respect to 
v): It consists of all permutations a € Sym(p) for which there exist matrices 
A G GL(C) such that Avi = ^^-(j) for 1 < i < p. Clearly, Lmy{C) is always 
isomorphic to a subgroup of GL(C), which we call GL^(C). In the particular case 
of a polytope P with an associated cone C generated by {v'^ = (1, Wj))i<t<p the 
group GL„(C) is actually the group of affine transformations preserving P. 

For every polyhedral cone C and every set of generators v we have 

Lin„(C) < Proj(C) < Comb(C). 

Both inclusions can be strict and we can define the corresponding notions of equiv- 
alence. In Figure 1 we give examples of those notions for the cube. 

In Section 4 we prove that the projective symmetry group Proj(C) can be real- 
ized as Lin„(C) for a suitable choice of vectors v. However, in [4, 12, 26] some 
polytopes whose combinatorial symmetries cannot be realized as projective sym- 
metries are given. 

Using the implementation in [34] we have compared the linear, projective and 
combinatorial symmetry group of 4313 polytopes available from the web page of 
A. Paffenholz [43]. For these examples, only in one case is the projective symmetry 
group larger than the linear symmetry group. This example was the one obtained 
by applying the construction E2 [26] to the 4-simplex; it is projectively equivalent 
to the dual of the Johnson polytope J(5,2). For 75 of the 4313 examples, the 
combinatorial symmetry group is larger than the projective symmetry group. The 
additional symmetries are in most cases a factor of 2 but reached in two cases a 
factor of 36. 
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3. Computing Lin^,(C) 



In this section we give algorithms to compute the hnear symmetry group Lmy{C) 
and certain subgroups occurring in applications like combinatorial optimization. 

The linear symmetry group is easier to compute in practice than either the pro- 
jective or combinatorial group, at least in the typical situation where we have only 
a generator representation (or only an inequality representation) for the input. Fur- 
thermore, the computation of the linear group is used as a subroutine in our algo- 
rithms to compute the projective and the combinatorial symmetry groups. 

One practical method for computing Lmy{C) is to define a positive definite ma- 
trix 



and the graph G{v) on p vertices {1, . . . ,p} with vertex and edge colors Wij = 
vfQ~^Vj. As shown in [5] the linear group Lin^(C) is then the group of per- 
mutations that preserve the colors. In the particular case of a poly tope P C M"^ 
generated by (fi)i<i<p with associated polyhedral cone C generated by {v'^ = 
(1, f j))i<i<p the matrix allows to define a Euclidean scalar product on M"^ 
for which GL^(C) is the group of affine isometrics. 

By using the methods presented in Section 3.4 one can compute Lin„(C) for 
polytopes with a few thousand vertices. For polytopes with a large vertex set some 
reduction may be necessary. One idea used in [22], for which the poly tope has 
about 10^ vertices, is to compute the stabilizer of a vertex. 

In high dimensions a key bottleneck is the computation of the inverse of the 
matrix Q. One approach to the problem is to compute the inverse using double 
precision. A tolerance number tol has to be chosen and values of {Q~^)ij which 
are within tol have to be grouped. One then computes the automorphism group of 
the colored graph for the grouped colors and checks if the graph automorphisms 
can be represented by matrices. If they cannot, then tol has to be decreased or 
double precision is not enough. 

3.1. GL„ (Z) symmetries. The approach described in the previous section to com- 
pute Lin„(C) as a graph automorphism works reasonably well in most cases. How- 
ever, for polyhedra with many generators this approach may not be applicable be- 
cause the complete graph considered is too large. 

For certain polyhedra, the linear symmetries are realized as GL„(Z) matrices 
(integer matrices of determinant ±1). A particular class of polytopes that we en- 
countered with this property are so called consecutive ones polytopes (see, for 
example, [25]). These arise as the convex hull of m x n matrices with 0, 1 en- 
tries satisfying a consecutive ones property. Because these polytopes have about 
2m-n vertices in dimension mn, the graph construction for computing symmetries 
is infeasible even for small m and n. The linear symmetries can nevertheless be 
obtained very quickly in small dimensions using the approach discussed in this 
section. 



p 



(1) 




i=l 
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We now describe a method to compute the group GL^,(C,Z) := GLy{C) n 
GL„(Z). As discussed above, in some cases the induced permutations gener- 
ate Lilly {C). Let the representatives (i'i)i<i<p be integral. We conclude that 
Gq = {A e GL„(Z) : A^QA = Q} (with Q as in (1)) contains GU(C,Z). 
In particular we obtain GLy{C, Z) as the setwise stabilizer of . . . , Wp} in Gq. 
The group Gq is the automorphism group of a lattice, that is of a discrete additive 
subgroup of the integral vectors. The matrix Q is an integral gram matrix of a 
lattice basis. The group Gq can thus be computed with the algorithm of Plesken 
and Souvignier [28]. It remains to perform a stabilizer computation in this matrix 
group. In principle, the Plesken-Souvignier algorithm can be adapted to include 
this stabilization. As the matrices A are generated row by row in a backtrack algo- 
rithm, we can check whether the current set of rows of A stabilizes a projection of 
{f 1, . . . , Vp} accordingly. If it violates this stabilization property, we may discard 
the entire candidate branch. 

The Plesken-Souvignier algorithm computes a set 5 C of short lattice vec- 
tors, which may be a difficult task for itself. Then a backtrack search on S is 
employed to compute Gq. If we briefly ignore the computational cost of S, the 
lattice approach has the advantage to work with an n x n instead ofapxp matrix. 
Thus for polyhedra which are generated by many rays and for which S is not too 
difficult to compute this may be a viable alternative. 

Theoretically, this method could be used to compute the linear symmetries of 
any cone C. In particular, if {vi, . . . , Vp} C and contains the standard basis 
ei, . . . ,en, the groups GLt,(C,Z) and GLy{C) are the same. If this condition is 
not satisfied by the input, we may transform C linearly without altering its linear 
symmetries. 

In some applications the goal is to find some GL„(Z) subgroup of GL„(C), 
rather than the full hnear symmetry group. One such application occurs in [10] 
when computing the Delaunay tessellations of a n-dimensional lattice L. We 
denote by Isom(P) the group of isometrics of a Delaunay polytope P and by 
Isomi(P) the group of isometrics of P that also preserve the lattice L. One sim- 
ple technique to find Isouil^P) is to iterate over Isom(P) and keep the elements 
belonging to lsomL{P). The finite group Aut(L) of isometrics of L preserving 
is identified with the group of isometrics of the quotient M"/ l. The center c of 
the empty sphere of P is expressed as ^ with f G L and the group Isomx, (P) is 
identified with the stabilizer of c G 1^"/-^^ by Aut(L). For any divisor d of m we 
can consider the stabilizer in K"/ L' with L' = L/d and the factorization of m gives 
a sequence of stabilizers that converges to Isomi(P). However, the most powerful 
technique to compute lsom.L{P) in practice is to apply the Plesken-Souvignier al- 
gorithm to the following homogenization: We define a (n + 1) -dimensional lattice 
L' spanned by all (0, v) and (0, v — c) for v G L. The automorphism group of 
L' contains an index 2 subgroup isomorphic to Isomi,(P), see [10] for more de- 
tails and [34] for an implementation of this technique which is very similar to the 
previous use of the Plesken-Souvignier algorithm. 
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In the general case of a polytope P for which we want the whole group GLt, (P, Z) 
the above techniques do not apply. Apart from the method of iterating over all ele- 
ments, one technique could be to add elements to the vertex set of P until we get a 
set that spans Z". 

3.2. Symmetries of Integer Linear Programming Problems. In integer linear 
programming, one optimizes over the intersection of a polyhedron with the integer 
lattice. From a mathematical point of view, the natural symmetries thus preserve 
the polyhedron and the integer lattice, i.e. are GL„(Z) symmetries. As far as we 
know, these general symmetries are not used in existing integer optimization soft- 
ware. Instead, the common practice is to consider only coordinate symmetries, 
i.e. permutations of coordinates that are automorphisms of the polyhedron. These 
symmetries turn out to be easier to compute, and also straightforward to work with 
in integer programming solvers. 

Several authors have been concerned with ways to compute coordinate symme- 
tries, by reducing the problem to a graph automorphism problem (see [29, 31,3, 
23]). Coordinate symmetries are isomorphic to the automorphisms of the following 
bipartite graph. Its vertex set is the union {vi, . . . , Vp} U {xi, . . . , Xn} of gener- 
ators (coming from inequalities) and variables. Between each pair Vi and xj we 
add an edge colored by the coefficient of variable Xj in generator Vi. The input is 
transformed to a bipartite edge colored graph, which is simpler than the complete 
colored graph required for more general symmetries. In integer programming we 
also have an objective function, which has to be considered for symmetry compu- 
tation. We can deal with this by coloring the graph vertices that correspond to the 
variables xi, . . . , x„ by their respective objective coefficient. 

Using this graph and transformation techniques detailed in Section 3.4, the per- 
mutation symmetries of 353 polytopes from mixed integer optimization were com- 
puted in [27]. Some of the larger instances, which had more than one million 
variables or facets, were still computationally tractable. In 208 polytopes a non- 
trivial symmetry group was found. For the 50 smallest problems, with dimension 
less than 1500, we computed also the linear symmetry group. We found that in 
6 out of these 50 cases the linear symmetry group is larger then the permutation 
symmetry group. All these linear symmetries are realized by integral matrices. 

3.3. Centralizer subgroups. We now give an algorithm for centralizer subgroups 
that is useful particularly in applications in the Geometry of Numbers; several 
examples follow the proof. For a given set B C Af„(]R) we want to find the group 

GU{C, B) = {Ae GL^(C) s.t. AB = BA for B £ B} . 

Without loss of generality we may assume that the set B is linearly independent, 
contains the identity matrix and is written as {Bi, . . . , Br} with Bi = In- We 
denote by Lmy{C, B) the corresponding isomorphic permutation group which is a 
subgroup of Lilly (C). 

Theorem 1. If B = . . . , i?,.} is a set of n x n-matrices with Bi = In then 
the group Lin^(C, B) is the group of permutations a preserving the directed colored 
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graph with edge and vertex colors 

"^ij = {vjQ^^BiVj, . . . ,vfQ^^BrVj) 

with 

p 

g = ^ Vivf 

i=l 

Proof. If A G G'Ly{C,B) then ABi = BiA and Avi = v^jiiy Hence one gets 
by summation that AQA^ = Q or equivalently Q^^ = A^Q~^A (obtained from 
Q^^ = {A'^)~^Q~'^A^^ by left and right multiphcation). So, if one writes h^j = 
vfQ~^BkVj then one gets 

h'^j = vfQ-^BkVj 

= vjA^Q-^ABkVj 
= {AvifQ-^Bk{Avj) 

So, any A induces a permutation of the p vertices preserving the vector edge 
color Wij. 

Suppose now that a S Sym(p) satisfies Wij = Wfj(^i-)cr{j)- Then, since Bi = In 
we have vfQ^^Vj = i'J(j)Q~^i'cr(j)- By an easy linear algebra computation (see 
[9, 5] for details) we get that there exists A G GL„(M) such that Avi = v„(^iy If 
one writes Wi = Q^^Avi then 

wfABkVj = vjA^Q-^ABkVj 
= vjQ~^BkVj 



a{i)a{j) 

= vjA^Q-^BkAvj 
= wfBkAvj 

Since the families (vj) and (wi) span M" we get AB^ = Bj.A. □ 
There are many contexts where the above theorem is useful. For example if one 
wishes to find the group of elements A € GL„(Q[-v/5]) then one way is to express 
the elements as elements of GL2n (Q) that commute with the multiplication by \/5. 
This is very useful when working with Humbert forms [1] whose symmetry group 
in GL„(Z[\/5]) correspond to a small subgroup of the full group in GL2n(Z)- 

Another such example is if one wishes to find the elements belonging to GL„ (H) 
with IHI the Hamilton's quaternions. GL„(]HI) acts on H" by multiplication on the 
left and it is characterized in GL^n (K) by the fact that it commutes with the Hamil- 
tonian multiplication on the right. 

Another example is if one wishes to compute the group GLy{C, W) of elements 
of GL„(]R) preserving a polytope C and a vector space W. Any such element will 
be an isometry for the scalar product defined by Q^^ in the proof of Theorem 1 
and so will commute with the orthogonal projection on W. 
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One would wish for a similar characterization for elements of GLv{C) that pre- 
serve a set B setwise by conjugation and so compute normalizer groups. But we 
do not think that a similar characterization is possible. 

3.4. Computing with vertex and edge colored graphs. Many of the algorithms 
presented above depend on the computation of the automorphism group of a graph 
whose vertices and/or edges are colored. The complexity of the graph isomorphism 
problem is uncertain as it is one of the rare problems in NP which is neither known 
to be NP-complete nor in P. 

For practical computation there exists graph isomorphism software (see [37, 36, 
33, 39]) that can compute automorphism groups of graphs. Such programs usually 
use the partition backtrack algorithm and can compute the automorphism groups 
of large graphs but their run time is exponential in the worst case. These programs 
usually cannot handle edge colored graphs and suffer from a performance penalty 
when using digraphs. Hence, one needs reduction techniques. 

There are several techniques for reducing an edge colored graph to a vertex 
colored graph (i.e. a complete graph with only two edge colors). One transfor- 
mation described in [31] replaces each c-colored edge {a, b} with an intermediate 
c-colored vertex m, which has edges connecting it to a and b. The obtained graph 
has n + n(n — l)/2 vertices which makes the transformation expensive. For the 
bipartite edge colored graphs that occur in Section 3.2 this can be improved by an 
idea given in [29] . Instead of adding intermediate vertices for all edges, we com- 
bine some of those with the same color. Define the bipartition as [S, S'). For each 
i £ S let c ^ S' be the set of vertices which are incident to i with an edge of 
color c. Then it is enough to introduce an intermediate c-colored vertex m with 
edges to i and to all elements of Xj c- For many integer optimization problems 
these sets c are often large, thus the number of vertices in the graph is usually 
substantially reduced. If l-Sj > |S"|, it may be advantageous to combine edges the 
other way around. 

For general edge-colored graphs we use the following method proposed in the 
user manual of nauty, [37]. Suppose that we have M colors. Then any color 
can be expressed as a 0/1 word of length [log2(M)]. Hence, the automorphism 
group can be obtained from the superposition of [log2(M)] vertex colored graphs 
and so from a graph with p[log2(Af )] vertices. This solution has good complexity 
estimates but the preceding method is often the best for the bipartite graphs from 
integer optimization (see [27]). 

Colored digraphs can be transformed into colored graphs with twice the number 
of vertices. A vertex a corresponds to a pair {a, a'} and a directed edge (a, b) to an 
undirected edge (a, b'). Let c and c' two colors that do not occur as directed edge 
color. We assign edges (a, b), respectively (a', b') the color c, respectively c'. 

4. Computing Proj(C) 

For a given polyhedral cone C, the group Proj(C) is the group of permutations 
of extreme rays that are induced by a linear transformation of C. We give a method 
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allowing the computation of the projective group in practice. It is based on a de- 
composition for polyhedral cones and on some Unear algebra tests. 

A polyhedral cone C generated by rays E = {Ri)i<i<p in a vector space V is 
said to be decomposable if there exist two non-empty subspaces Vi, V2 such that 
E = {EnVi)\JiEn V2) and V = Vi® V2. 

Theorem 2. For a polyhedral cone C generated hy E = {Ri)i<i<p in a vector 
space V, there is a unique decomposition 

V = ®i<k<hVk 

with 

E = Ui<k<hE n Vk 

and the cone Ck generated by EnVk being non-decomposable. This decomposition 
can be computed in time O (pn^) . 

Proof. The existence of the decomposition is obvious since the vector space V is 
finite dimensional and so the decomposition process has to end at some point. The 
uniqueness follows immediately from the fact that if there are two decompositions 
by Vk and V/ then the family of subspaces Vj. n also defines a decomposition. 

The method for computing a decomposition is the following. First select some 
generators Vi of Ri. Then find (e.g. by Gaussian elimination) a n-element set 
S C {1,.. . ,p} such that {vi)i^s is a basis of V. Every vector Vi fori < i < pis 
then written as Vi = J2keS '^k,iVk- The supports 

Si = {keS : ak,r ^ 0} 

determine the edges of an hypergraph on n vertices. The connected components 
of this hypergraph correspond to the summands of the decomposition of V. Both 
computing S and finding the connected components can be done in O (pn) time. 

□ 

The above decomposition allows one to determine the projective symmetry group. 
We say that two cones have the same the projective isomorphism type if there is a 
bijective linear map between them. 

Theorem 3. For a polyhedral cone C generated by rays E = {Ri)i<i<p in a vector 
space V, decomposed into V = ®i<k<hVk< define Ck to be the cone generated by 

EnVk. 

If the projective isomorphism type j tor 1 < j < q occurs nj times among the 
Ck then we have the equality 

Proj(C) = Wr(Sym(n,),Proj(Cfc,)) 

where Wr stands for wreath product and Ck^ a representative for the jth type. 

Proof. Suppose that we have an element / G Proj(C). This element permutes the 
Ck but must also preserve the projective isomorphism types. Hence, it belongs to 
the above mentioned product. The reverse inclusion is trivial. □ 
We now expose a method for computing the projective symmetry group of a 
non-decomposable polyhedral cone C. Combined with the above theorem, this 



SYMMETRY GROUPS OF POLYHEDRA 



11 



will allow US to compute the projective symmetry group of arbitrary polyhedra. 
We first present the following structural result: 

Theorem 4. Suppose C is a non-decomposable polyhedral cone generated by rays 
iRi)i<i<p- Then 

(i) We have the isomorphism 

GL(C) ^R^x Proj(C). 

(ii) There exist vectors Vi such that Ri = M+Vj and 

Lin^(C) = Proj(C). 

Proof. To show (i), we have to determine the kernel K of the homomorphism 
GL(C) ^ Proj(C) and show that K splits from GL(C), i.e. GL(C) ^ K x 
Proj(C). For A ^ K we have ARi = aiRi with > 0. If the values ai were not 
all the same then this would automatically give a decomposition of the space (since 
each class of rays with the same multiplier will be subdimensional) contradicting 
the non-decomposability of C. Thus K = {A/ | A G M^}. Therefore GL(C) is a 
central extension K. Proj(C) of Proj(C) by K, i.e. K is normal in GL(C) and in 
particular lies in the center of GL(C) (see e.g. [13, (11.8)]). A classical theorem 
due to I. Schur states that every non-split central extension (i.e. a central extension 
A.B which is not isomorphic to ^4 x i3) of a finite group G can be obtained as a 
homomorphic image of M{G).G, where M{G) is a finite Abelian group, called 
the Schur multiplier (cf. [30] or [13, (11.17)]) of G. Thus K > M', where M' is 
a quotient group of the Schur multiplier M(Proj(C)) of Proj(C). As K is torsion- 
free (i.e. has no non-identity elements of finite order), this implies that M' = 1, 
and GL{C) = K X Proj(C), as claimed. 

To show (ii) we can, by (i), identify Proj(C) with a subgroup H of GL(C). This 
identification is unique, as Proj(C) is the subgroup of elements of finite order. For 
each Proj(C) -orbit Ok of extreme rays we choose a representative ray Rk € Ok 
and a vector Vk such that Rk = R+Vk- For each element R £ Ok we can find an 
f £ H such that R = f{Rk). If there is another f £ H such that R = f'{Rk) 
then h{vk) = Cvk with h = f^^f and C > 0. Since /i G if, it has finite order, 
and thus G = I. This means that f{vk) is uniquely defined. It is then clear that for 
this choice of generators Proj(C) = Lin„(C). □ 

Remark 5. The conclusions of Theorem 4 cease to hold in a more general setting 
of a configuration of points in a projective space over C. In the following we 
construct a counterexample. 

First we take the group G characterized as 2 ■ (namely, number 28 in GAP 
database of small groups [35]) and its faithful 2-dimensional representation over C. 
The center of G in this representation is it Id2. We then take an element of order 8 
in G and compute one of its eigenspaces, corresponding to an 8-th primitive root of 
unity. There are 6 images of the eigenspace under G. Thus we obtain a transitive 
action of G on a 6-tuple of lines and so on 6 points in P^{C). The action on the 6 
lines defines a group 5*4. But G is not isomorphic to 2 x 5*4. 
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Theorem 6. Suppose C is a non-decomposable cone generated by rays {Ri)i<i<p 
in M". Testing if a permutation a G Sym(p) belongs to Proj(C) can be done by 
solving a linear system with np equations and p unknowns. 

Proof. Let us slightly abuse notation and denote by = (^'^(i), • • • the 
matrix with columns being a set of generators for the Ri, i.e. Ri = M_|_t!j, for 
1 <i < p, permuted by a permutation ji. Respectively, let Ufj, denote the submatrix 
of Vn consisting of its first n columns, and V := Vi^, U := Ui^. Without loss of 
generality, U is invertible. 

The sought matrix A G GL(C) — a preimage of a — must satisfy Avi = aiV^(^i^, 
ai > 0, for each 1 < z < p. In the matrix form this can be written as AV = 
Va diag(ai, . . . , a^). In particular, AU = 11^ diag(ai, . . . , q;„), implying 

(2) A = Uadi&.g{ai,...,an)U-'^. 
This implies 

(3) [/(^diag(ai, . . . ,Q;„)i7~V = Vad\a,g{ai, . . . ,ap). 

This is a homogeneous linear system having np equations and unknowns ccj, for 
1 < i < p. Denote by SV the solution space of (3). A solution a is acceptable 
if and only if a > (the latter implies det(A) 7^ by (2)). These conditions are 
open conditions, so if there is one such solution then there is an open ball of such 
solutions of dimension q := dimSV, as well. But we know by Theorem 4 that 
g < 1 for non-decomposable cones. Thus either q = 1, or a ^ Proj(C). If gr = 1, 
we can find a nonzero solution a of (3) and test that ±a > 0. If there is no such a, 
we conclude that a ^ Proj(C). Otherwise, picking the right sign of a, we find A 
using (2) and conclude that a € Proj(C). □ 

Theorem 6 gives a constructive way to compute Proj(C). Combining with the 
intermediate subgroup algorithm to compute Comb(C) given in Section 5 gives a 
more practical method to compute Proj(C). An easy situation is when Lini,(C) = 
Comb(C), which of course implies Lin„(C) = Proj(C). 

Let us take ai > 0. The group LinQ,^(C) depends on a and is a subgroup of 
Proj(C). By Theorem 4 there exist a such that LinQ,„(C) = Proj(C) and it is 
interesting to know when equality occurs. For any a > the group LinQ,^(C) is 
a symmetric group acting on p points which defines an orbit partition OP (a) of 
{!,..., p}. 

Theorem 7. (i) If {ai)i<i<p and {a^)i<i<p are two sets of positive multipliers and 

OP(a) = OP(a') then Lin^^C) = Lin«/^(C) 

(ii) If (ai)i<j<p is a set of positive multiplier and LinQ;^(C) is transitive on 
{l,...p}then Lin«^,(C) = Proj(C). 

Proof. Let us prove (i). We decompose C into non-decomposable components 
for 1 < k < h and denote by the corresponding subset of {1,. . . ,p}. Let us 
take an orbit O under Ga = Lino,^(C). If x e O D Sk and H = Stabc^ (Sk) then 
Theorem 3 giving the expression of the projective symmetry group in terms of a 
wreath product is also valid for the linear automorphism group. This implies that 
the orbit of x under H is exactly OCiS}-. Our assumption 0P(q:) = OP(a') imphes 
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that O n Sk is also an orbit under H' = Stabc^, (5^). Furthermore, by the non- 
decomposabiUty of Ck we know that ai is determined up to some constant factor 
for Z G On Sk- Thus there exists a /3 > such that ai = jia'i for / G O n 5^. This 
implies that the linear automorphism groups of for {aiVi)i^Sk ^^^^ {a[vi)i^Sk 
are equal. Since OP(a) = OP(a') we know that the isomorphism type of the 
component Ck under Ga and G^i are the same. Since both groups are actually 
direct products of wreath products they are necessarily equal. 

To prove (ii), note that by Theorem 4 there exist some multiplier a' such that 
Ga' = Proj(C). Ga is a subgroup of Ga' so OP(a) is a partition induced from 
OP (a') by splitting some orbit. But OP (a) is reduced to only one component so 
OP(a) = OP(q') and one concludes. □ 

By using item (ii) above one can conclude in some cases that the linear group is 
actually the projective group. 

5. Computing Comb(C) 

Recall that Comb(C) is the maximal symmetry group of a polyhedral cone that 
preserves the face lattice. For many polyhedral computations, this is the largest 
group of symmetries that can be exploited. Although no efficient methods are 
known for the general case, in this section we describe some techniques that can be 
useful in certain practical computations. The general idea is to construct a "sand- 
wich" Gi < Comb(C) < G2 between groups Gi and G2 that are easier to com- 
pute. We present a technique based on double coset decomposition that can be used 
to speed up the testing of such a group inclusion for strictness. 

We define the group of combinatorial symmetries Skelfc(C) to be the group of 
symmetries preserving the faces of dimension at most k. In particular Skeli(C) = 
Sym(p), Skelfe+i(C) < Skelfc(C) and Skeln-i(C) = Comb(C). 

For a chosen set of generators (fi)i<i<p and an integer /c > we have the 
inclusion 

(4) Lin^(C) < Proj(C) < Comb(C) < Skelfc(C). 

Assuming that we know the set Tk of /c-dimensional faces, the group Skelfc(C) 
is isomorphic to the automorphism group of a vertex colored graph on p + \J^k\ 
vertices. The reason is that if an automorphism preserves all the A;-dimensional 
faces, then it preserve all the intersections and so all the faces of dimension at most 
k. The A;-dimensional faces F are thus described by the set S of vertices contained 
in them and so we can build a bipartite graph onp+\Tk\ vertices that encodes this 
relation. 

By the chain of inclusions in (4), the group Comb(C) is located between two 
groups which are both automorphism groups of colored graphs. If we can prove 
that for some ko we have Lm.u{C) = Skel^p (C) then we conclude that Comb(C) = 
Lm^{C) and we are finished. This is the most common method (see [8, 7]) for 
computing combinatorial symmetry groups: determine the set of faces of dimen- 
sion at most ko, determine Skel^g (C), test if the elements of Skelfco (C) are actually 
in Lin^, (C) and if yes obtain Skelfc,, (C) = Comb(C). But this does not always work 
since in some cases Lin„(C) / Comb(C). One case where it is guaranteed to work 
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is for simple polytopes, i.e. ones for which every vertex is adjacent to exactly n 
other vertices, for these polytopes Comb(C) = Skel2(C) [2, 15, 11]. 

Typically, the number of A;-dimensional faces becomes impractically large for 
some intermediate values of k. An alternative method for computing Comb(C) 
is simply to compute the whole set of facets and then compute Skel„_i(C) = 
Comb(C) directly. The problem of this method is that the set of facets may be too 
large for this approach to work and sometimes the facets precisely what we want 
to compute in the end. 

Suppose now that we have a groups Gi, G2 with Gi < Comb(C) < G2, We 
now consider a method based on double-coset decomposition to test if these inclu- 
sions are strict, possibly replacing Gi with a larger subgroup of Comb(C) in the 
process. 

Let us assume that we have computed the orbits of facets of C up to Gi. The 
possible methods for doing such a computation are reviewed in [5]. Denote by Oi, 
. . . , Or the orbits of facets, for which we select some representative Fi, ... , F^. 
They are encoded by their vertex incidence as subsets 6*1, . . . , 5"^ C {1, . . . 

The first step is to be able to say whether a permutation a of the vertex set of 
C does belong to Comb(C). A permutation a G Sym(p) belongs to Comb(C) 
if and only if any image a{Si) is in a Gi-orbit of one of our representatives Sj. 
Such in-orbit tests are done using permutation backtrack algorithms [20, 21] that 
are implemented, for example in GAP [35] and PermLib [39]. 

Now suppose we are given three groups Gi C H C G2 with H described by 
an oracle (e.g. the procedure based on Gi-orbits just described). We propose an 
intermediate subgroup algorithm for determining an explicit representation for H. 
We can do a double coset decomposition of G2 using the subgroup Gi : 

s 
i=l 

with Qi G G2 and GigiGi n Gigi'Gi 7^ if and only if i = i'. Suppose that 
g ^ H, then since Gi C G2 for every /, /' G Gi, we have fgf G H. So, for a 
given g G G2 either GigGi C H or GigGi r\ H = %. This allows a reduction in 
the number of oracle calls. Additionally if we found a g £ G2 — Gi that belongs 
to H then we can replace Gi by the group generated by g and Gi and recompute 
the double coset decomposition. In particular this method allows us to compute 
Comb(C) without having to iterate over all elements of Skelfc,, (C) and test whether 
they belong to Comb(C). 

It is not clear how to do much better since in general one needs the facets in 
order to get Comb(C). The underlying assumption to get good performance using 
the intermediate subgroup algorithm is that the index [G2 : Gi] is not "too large", 
i.e. [Skelfc(,(C) : Lin^,(C)] is not "too large". 
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